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ABSTRACT 

We study relations between the X-ray luminosity, orbital period and absolute near- 
infrared magnitude of persistent low-mass X-ray binaries (LMXBs) . We show that of- 
ten optical and near-infrared spectral energy distribution of LMXBs can be adequately 
described by a simple model of an accretion disc and a secondary star reprocessing 
X-ray emission of a central compact object. This gives us an evidence that using an 
X-ray luminosity and an absolute infrared magnitude of a persistent LMXB one can 
make reliable estimate of its orbital period. Using a sample of well-known LMXBs, we 
have constructed a correlation of Lx, Porb and Mjc values which can be approximated 
by a straight line with the RMS scatter at the level of ~0.3 mag. Such a correlation, 
being to some extent an analogous to the correlation, found by |van Pa radijs fc Mc-| 
Clintock (1994), might be helpful for future population studies especially in the light 
of forthcoming surveys of the Galaxy in X-ray and infrared spectral domains. 
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1 INTRODUCTION 

Low-mass X-ray binaries are binary systems with compact 
objects (neutron stars or black holes) accreting matter from 
a low-mass secondary companion. First such systems were 



discovered about 50 years ago ( Giacconi et al. 1962 ) and 



presently more than a hundred of them are known in our 
Galaxy (e.g. Liu, van Paradijs, & van den Heuvel|2007[> and 



Fabbiano 



hundreds of them in external galaxies (see e.g. 
|2006| |Evans et al.|2010| . Population of low- mass X-ray bina- 
ries provides us with tools to study different physical mech- 
anisms, that define the population features (e.g. the energy 
transfer efficiency during the common envelope stage of the 
binary, the gravitational wave emission influence, the aver- 
age magnetic field strength in old neutron stars), and also 
some properties of a host galaxy (e.g. the age of its stellar 
population). In order to reach these goals one needs: (1) to 
better understand the laws of formation and evolution of 
LMXBs, and (2) to gather a clear sample of LMXBs with 
known parameters (orbital periods, X-ray luminosities, etc). 

One of the major obstacles with the latter item is the 
lack of measurements of orbital periods for the majority of 
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known LMXBs. This is mainly caused by the fact that the 
orbital modulation of the X-ray emission is rarely observed 
due to geometrical reasons (small companion star covers in- 
significant fraction of the sky for the compact object, there- 
fore very few known LMXBs demonstrate eclipses and/or 
dips), whereas the detection of the optical or infrared (IR) 
modulation of LMXBs emission is challenging due to ex- 
treme faintness of these objects. 

However, it is sometimes possible to estimate the LMXB 
orbital periods indirectly, making use of the fact that their 
optical and IR emission originates mainly from reprocess- 
ing of the central X-ray emission by the outer accretion disc 
(see e.g. McClintock et al.|1979 van Paradijs & McClintock 
1994|) . The efficiency of this approach was demonstrated by 

19941, who showed that the 



van Paradijs & McClintock 



absolute optical magnitude of LMXBs had very clear cor- 
relation with their X-ray luminosity and orbital period. It 
should be kept in mind, however, that this method is phys- 
ically justified for persistent LMXBs, which possess station- 
ary accretion flow with a constant incident X-ray flux. A 
transient-like behaviour might lead to the large variations 
of X-ray flux or/and the accretion flow structure and thus 
to the deviations from any possible correlation. One might 
also expect significant deviations if the central X-ray emis- 
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sion of the source is hidden from an observer due to the 
very large inclination angle of the binary system (so called 
Accretion Disc Corona/ADC sources like 4U1822-37). 

Unfortunately, the majority of LMXBs in the Galaxy 
(either transient or persistent) reside in the Galactic plane 
and thus suffer from severe dust extinction reaching Ay ~ 
10 — 30 mag making them completely unobservable in the op- 
tical domain, thus shifting observational methods to longer 
wavelengths. Due to the fact that reliable measurements of 
LMXBs absolute magnitude in the infrared spectral bands 
can be done for much more sources, it is then reasonable to 
calibrate diagnostics similar to that of van Paradijs & Mc- 
|Clintock| |l994) in the IR, in particular, in the K spectral 
band. 

However, spectral energy distribution (SED) of sources 
at these wavelengths might be contaminated by optically 
thin synchrotron emission, sometimes observed in the far 



end of NIR up to radio range in persistent ( jMigliari et al 



20101 and especially often in transient LMXBs (see e.g. 



Cor- 



bel & Fendei 


2002) |Russell et al. 


2006| |Russell, Fender, & 


Jonker 


2007 


Shahbaz et al. 


2008 


Russell & Fender||2008|. 



This issue is to be carefully explored before we are able to 
construct reliable period-luminosity relation for the K band. 

In this paper we show that the general shape of SED 
of persistent LMXBs in the optical-NIR spectral range can 
be reasonably well described in the framework of a model, 
where the optical-NIR emission comes from the reprocess- 
ing of the central X-ray flux in the accretion disc and the 
secondary star. This gives us an evidence that there is a 
physical connection between the main properties of the bi- 
nary and its infrared luminosity. We present the correlation 
of an absolute magnitude of LMXBs with their X-ray lumi- 
nosity and orbital period similar to the one found in the V 



band by van Paradijs & McClintock ( 1994 1. In this work we 
consider only persistent LMXBs to be sure that their accre- 
tion disc structure does not change much at different epochs. 
We consider only LMXBs which harbour neutron stars as a 
compact object in order to reduce the uncertainties related 
with the mass of the compact object. 



2 LMXBS SPECTRAL ENERGY 

DISTRIBUTION IN THE OPTICAL AND 
NIR 

In order to check the validity of our approach to interpre- 
tation of the LMXB luminosity in optical and NIR spectral 
domains, we have studied broad band spectral energy distri- 
bution of the three binary systems, Cyg X-2, Sco X-l, and 
4U0614+091, which can be considered as illustrative. These 
three systems cover wide range of orbital periods, from 0.81 
to 236.27 h, and have more or less reliable distance, inter- 
stellar extinction and in two cases (Cyg X-2 and Sco X-l) in- 
clination estimates. Their parameters, adopted in this study 
for spectral modelling, are presented in Table [l] 

It is important to keep in mind that in some cases a con- 
siderable contribution from an optically thin emission can be 
seen at different states of X-ray binaries (see e.g. |Russell et] 
al.|2006 1. In particular, detection of polarized infrared emis- 
sion of sonic X-ray binaries, exceeding levels expected from 
interstellar scattering, hints at synchrotron origin of some 
fraction of their infrared flux (see e.g. Shahbaz et al.|2008 |. 



Therefore, in order to test the applicability of our model 
that accounts for the optically thick emission, we have inten- 
tionally selected X-ray binaries at different spectral states, 
namely high/soft (Sco X-l and Cyg X-2) and low/hard 
(4U0614+091) state. 

For calculation of spectral energy distribution of 



LMXBs we adopt a simple model (see e.g. |Tjemkes, van 
Paradijs, fc Zuiderwijk||1986| |Q'Brien et al.||2002| |Mesc~ 



eryakov, Revnivtsev, fc Filippova|2011 1 when all surface el- 



ements in a binary system (which comprises an accretion 
disc and a secondary star) emit black body radiation with 
the local temperature defined by X-ray irradiation and in- 
ternal heating. Shape of the star was assumed to be of Roche 
lobe geometry. 

Temperature of the illuminated side of the star was cal- 
culated taking into account heating due to incident X-ray 
flux: 

m4 _ rpi . ??*£xcos6>* 

where L x - luminosity of the central X-ray source, 8+ - angle 
between X-ray source direction and surface normal, d - dis- 
tance from X-ray source to the surface element, a - Stefan- 
Boltzmann constant, rji, - fraction of reprocessed X-ray emis- 
sion. Temperature of the part of the star not illuminated by a 
compact ob j ect wa s taken following |de JongTv an Paradijs, & 
Augusteijn| |l996| as T*, = 58OO(i/L ) 1/4 (i?/7? )- 1/; ' K 



and mass-luminosity and mass-radius relation from |Tout et| 
|al"| |l996). For the fraction of X-ray emission reprocessed 
and reradiated from the surface of the secondary star we 
assume 77* = 0.6 (see e.g. London, McCray, fc Auer||1981l. 



For the effective temperature of the disc surface element 
we used the common relation, which includes viscous heat- 
ing and heating by X-ray irradiation in the outer parts of 
geometrically thin accretion disc: 

3GA/!i x , 77 d L x 



Ti 



+ 



\ R )ou 



(n - 



(2) 



8nec 2 aR 3 ' 4naR 2 

We use Mi — 1.4 Mq and e = 0.1, - mass of the compact 
object and accretion efficiency, respectively, for NS LMXB. 

The fraction of the flux, intercepted by the accretion 
disc, was calculated assuming a shape of the intercepting 
surface as disk with height being a function of its radius: 
H (x R n , we assume n = 9/7 ( |Vrtilek et al.|1990| >. We note 
here that intercepting surface is not necessarily the optically 
thick accretion disc itself, it might be the corona that inter- 
cepts the X-ray flux and redirects it to the optically thick 
accreti on disc, see e.g.[jTmenez-Garate, Raymond, fc Liedahl| 
( |2002[ ); pescheryakov, Shakura, fc Suleimanov| ( |2011[ ). We 
assume ratio of hight of the intercepting surface to its radius 
at the outer edge of the disc (H/R) ou t = 0.1. Parameter r/d 
- the effective fraction of emission, reprocessed in the ac- 
cretion disc-corona system (ratio of X-ray flux thermalised 
in optically thick disc to the incident X-ray flux) is not yet 
accurately known, different authors estimate its value from 
less tha n 0.1 (|de Jong, van Paradijs, fc~A ugusteij n|1996| l up 
to ~0.5 ( Vrtilek et al.|1991 1. In our work we adopt the value 
rj d = 0.25. 

The accretion disc radius around the compact object 
was taken to be equal to the tidal radius as estimated by 



Paczynski (19771. More detailed description of our model 
can be found in [Meschcryakov, Revnivtse v, fc Filippova] 
(20111. We would like to note here that our ability to pre- 
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Table 1. Parameters of the three binary systems, for which we 
have calculated optical-NIR SEDs. Orbital periods were adopted 
from |Ritter fc Kolb| p003} . Cyg X-2 and Sco X-l, being members 
of the so called Z-sources, were assumed to have the Eddington 
luminosity, whereas the luminosity of 4U0614+091 is taken from 
|Revnivtsev et aL| l |201l} . See the text for other parameters, com- 
mon for these systems. 



Table 2. Photometric measurements of Cyg X-2. We have as- 
sumed its interstellar extinction to be Ay = 1.34. 





Cyg X-2 


Sco X-l 


4U0614+091 


Distance(kpc) 


11.6 iO.3 1 


2.8±0.3 2 


3.2 3 


-^orb 


236.27 


18.94 


0.81 


Mi 


1.71 4 


1.4 


1.4 


1 


0.34 4 


0.3 s 


0.04 6 


log L x (erg/sec) 


38.3 


38.3 


36.5 


Inclination 


62. 5 7 


38 5 


60 s 


A v 


1.34 9 


0.91*0 


2 11 



(1) -|Smale|l|1998[l, (2) -|Bradshaw, Fomalont, fc Geldzahler| 
l |1999|l, (3) -|Kuulkers et el!.| (|2010| ) , ( 4) -|Casares et a.I.| (|2010| ) , 
( 5) -|Steeghs fc Casares|l|2002| l, (6) - |Werner et al.| ( |2006| l, (7) 
— |Orosz fc Kuulkers| l |1999||, (8) - inclina tion is unknown, fixed 
at given val ue, (9) -|McClintock et al.| l|1984[ |, (10) - |Vrtilek| 



|et al.| l |199l) , (11) - |Migliari et al.| l |2010 l 



diet the absolute magnitude of LMXBs with the accuracy 
better than 0.5-0.8 mag (i.e. within a factor of two in the 
flux) mostly depends on the accuracy of our knowledge of 
the binary system inclination and the fraction of reprocessed 
X-ray emission (for fixed disc radius and shape of the X-rays 
intercepting surface). 



3 DATA 

In this work we consider only persistent LMXBs, which 
means that their optical, IR and X-ray fluxes do not vary by 
large factors. This is essential for present work, because we 
use here measurements of brightness of sources in optical, 
IR and X-ray energy bands, obtained not simultaneously. 

Distribution of X-ray fluxes of bright (with fluxes more 
than ~50 mCrab, i.e. which have enough statistics to be 
detected by All-Sky Monitor (ASM) of RXTE observatory 
in all single measurements) persistent LMXBs were anal- 
ysed by Postnov & Kuranov ( 2005 1 . It was shown there 



that the typical long-term root-mean-square flux deviations 
do not exceed ~20-30 per cent of their mean values. Long- 
term variations of optical emission of persistent LMXBs were 
studied in a number of works, in particular for sources Sco 
X-l ( |McNamara et al.|2005|>, Cyg X-2 (jO'Brien et al.|2004[ ), 



and 4U0614+091 ( |Hakala, Charles, fc Muhli||2011[ ) consid- 
ered in our paper. It was demonstrated that the long-term 
variations of their apparent optical magnitudes do not typi- 
cally exceed ~0. 1-0.3 mag (i.e. ~10-30 per cent flux varia- 
tion). 

Based on these evidences we assume additional ~10- 
30 per cent systematic flux uncertainty due to non- 
simultaneous nature of the measurements we use while 
analysing X-ray and optical luminosity correlations. 

To construct a broad band spectral energy distributions 
of our test objects we have used their multi-colour pho- 
tometry from the literature and also new publicly available 
data from Spitzer (instruments IRAC ( Fazio et al.|2004 | and 
MIPS ( [Rieke et al.|[2004| ) and WISE ( |Wright et aLpOlO] ) 
space observatories. 



Filter 


Vega mag 


A x 


Vega mag corr. 


U 


15.0 1 


2.05 


12.9 


B 


15. 3 1 


1.77 


13.5 


V 


14.8 1 


1.34 


13.5 


H 


14.0 2 


1.00 


13.0 


J 


13.39 3 


0.38 


13.01 


H 


13. 15 3 


0.24 


12.92 


K 


13.05 3 


0.15 


12.90 



Spitzer 



A, fim Flux, mjy 
4.5 1.4 
8.0 0.8 
24 8.5 



(1) - |Orosz fc Kuulkers|(l999"l l, (2) - USNO Bl, (3) - 2MASS 

Mid-infrared photometric measurements were obtained 
from the analysis of Spitzer basic calibrated data images. 
For each source we have determined the flux by summing 
counts in circular aperture with the radius of 6 arcsec (for 
Spitzer/MIPS at 24 /im we have used 9 arcsec aperture due 
to larger FWHM of the instrument PSF at these wave- 
lengths) centered at the object and subtracting the back- 
ground counts, collected in an appropriate field nearby. 

Another set of photometric measurements was ex- 
tracted from preliminary data release of the WISE satellite 
| |Wright et al.| [2010). Conversion of magnitudes into physi- 
cal fluxes was achieved by using the coefficients, presented 



Cutri et al. (2011). 



All photometric measurements were corrected for the 



interstellar extinction assuming the Rieke & Lebofsky ( 1985 1 
reddening law. 

Data at these long wavelengths are important for addi- 
tional check for the presence of optically thin synchrotron 
emission, sometimes seen in SEDs of LMXBs (see e.g. 
Migliari et al.|2010| >. 



3.1 Cyg X-2 

Photometric measurements of Cyg X-2 used in this work are 
given in Table [2] Model of the broadband SED, constructed 
assuming the above mentioned parameters is shown by solid 
curve. The model is not fitted to the data, but presented as it 
is, making use of the parameters, described in the text and in 
Table[l] Accurate position of the modelled SED with respect 
to Y-axis depends on the source distance, the fraction of 
emission reprocessed in thermal disc and the binary system 
inclination. 

It is clear from Fig. [I] that the adopted model ade- 
quately describes the shape of the broad band spectral en- 
ergy distribution of the source, except for its long wavelength 
range (8 and 24 /j,m). We make the following conclusions 
from this comparison: 1) the model well describes the SED 
of Cyg X-2 up to ~4 (im, and 2) at longer wavelengths we 
observe the indication of some additional emission compo- 
nent which is likely the optically thin synchrotron emission 




Figure 1. Extinction-corrected spectral energy distribution of 
Cyg X-2 in the optical and IR with a simple model of its emission 
(see text). The model is not fitted to the data, but presented as 
it is, making use of the parameters, described in the text and 
in Tab. [I] Dashed line represents contribution of the secondary 
star. Two sets of curves represent SEDs during maximum and 
minimum contribution of secondary star to the total brightness 
of the binary. Due to possible time variations of the optical/NIR 
emission of the source we have assumed 0.2 mag uncertainty of 
all photometric measurements. Open circles denote measurements 
obtained using Spitzer data. 



of non-thermal electrons. Polarization measurements, indi- 
cating that at these wavelengths the source might have a 
contribution from synchrotron component were presented 
by |Shahbaz et al.| ( [2008] ). 



3.2 Sco X-l 

A set of Sco X-l photometric measurements in the opti- 
cal and infrared, obtained at different epochs with differ- 
ent instruments, is presented in Table [3] It is clearly seen 
that Sco X-l exhibits significant variations of its spectral 
energy distribution at long wavelength end, while in the op- 
tical the changes are inessential (see e.g. |McNamara et alT] 



20051. It is not surprising because the source is known to 



be strongly variable at radio wavelengths (e.g. Pandey et 
|al.||2007| ), which means that it sometimes displays a spec- 
tral component emerging from non-thermal population of 
electrons. In addition to that it was shown that in the NIR 
Sco X-l sometimes exhibits polarization, exceeding the one 



induced by the interstellar dust scattering (Shahbaz etal 
2008| |Russell fc Fender|2008[ ), which can also be attributed 
to the emission caused by non-thermal electrons. Finally, 
yet another spectral component (different from the optically 
thick emission, considered by our model) can be found in the 
presence of hard X-ray tails, often observed from Sco X-l 
(e.g. |D'Amico et al.||2001| |Paizis et aL]|2006| l, and likely to 



Figure 2. The same as Fig. [T] but for Sco X-l. Open circles 
denote measurements obtained using Spitzer data, open squares 
- measurements from WISE data. 



be related to the emission of energetic non-thermal electrons 
in the accretion flow (see e.g. Migliari et al.|2007 |. 

However, in spite of these complications we note that 



one set of NIR measurements (namely those of Willis et al. 
1980 1 does show the Rayleigh-Jeans type SED in NIR sug- 



gesting that the contribution of an additional component at- 
tributable to non-thermal electrons at that epoch was small 
or negligible. We will therefore use this K magnitude esti- 
mate in our subsequent work. Importantly, the J — K colour 
can be used as an indicator of a presumably synchrotron 
component in the SED. In particular, our model predicts 
J — K ~ —0.15 mag for the source (see below), while the 
brightest NIR measurements reach J — K = 0.6 mag. The 
faintest NIR estimates give J — K — —0.19 mag, which is 
reasonably compatible with the predictions of our model. 



3.3 4U0614+091 



adopted 
from 



photometric 



We have 
4U0614+091 

seen (and it was mentioned in Migliari et al 
the 



Migliari et al. (20101 



measurements of 
It is clearly 
that 
at 



2010) 



source demonstrates optically thin component 
mid-infrared wavelengths and its contribution to K band 
can be comparable to that of the optically thick mechanism, 
considered by our model. For that reason we have calculated 
K magnitude of 4U0614+091 from J band measurement 
by Migliari et al. (|2010 1, using the model of optically thick 
emission (see Fig. |3|, and used this value in the subsequent 
work. 
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Table 3. Photometric measurements of Sco X-l. We have as- 



interstellar extinction to be Ay 
Filter Vega mag A\ Vej 


= 0.91. 

;a mag corr. 


B 


12.38 1 


1.78 


11.04 


V 


12.40 2 


0.91 


11.49 


R 


12. 3 3 


0.68 


11.62 


J 


11. 90 4 


0.25 


11.65 


H 


11. 54 4 


0.15 


11.39 


K 


11. 15 4 


0.10 


11.05 


J 


11.94 s 


0.25 


11.69 


H 


11.99 s 


0.15 


11.84 


K 


11.98 s 


0.10 


11.88 






Spitzer 





A, fim 

3.5 

4.5 

5.8 

8.0 

24 



Flux, mjy 
9.3 
7.1 
4.5 
.3.7 
1.3 

WISE 



3.35 
4.60 
11.5 
22.0 



10.4 
7.0 
2.5 
<3.0 



(1) - |McNamara et al.|j2005) , (2) - |McNamara et al. 1(2003) , 
(3) - USNO Bl, (4) - 2MASS, (5) - | Willis et al.| | |1980| l 




10" 10 5 
Wavelength, A 

Figure 3. The same as Fig. [I] but for 4U0614+091. 



4 Mk-S k RELATION OF PERSISTENT 
LMXBS 

4.1 Model 

Broad band spectral energy distributions, shown in the pre- 
vious section, provide us support that we adequately under- 
stand physical processes responsible for the observational 



appearance of LMXBs in the optical and NIR spectral do- 
mains. 

According to this model the absolute optical and NIR 
magnitude of LMXBs depend both on X-ray luminosity of 
the central source and on the size of the binary system, [van] 



Paradijs & McClintock (19941 showed that the luminosity 



of the system in optical V band scales approximately as a 
square root of its X-ray luminosity and as a power of 2/3 of 
its orbital period. Note that these indexes are theoretically 
expected for V band photometry if temperatures of matter 
that gives the main contribution to optical emission lie in 
particular range of temperatures (T « 10000-30000 K). For 
these temperatures the surface brightness of a black body 
emitter scales as Sv oc T~ 2 . For our case of the NIR K band 
the scalings should be different. 

In persistent LMXBs we expect that effective tempera- 
ture of accretion disc is everywhere above hydrogen recom- 
bination limit. The effective temperature at which a hot 
disc becomes thermally unstable is, according to |Dubus,| 
Hameury, & Lasota (20011, given by 

Here a is the standard Shakura-Sunyaev viscosity pa- 
rameter. 

For the temperature range T > Th the K spectral band 
(A ~ 2.2 fim) lies almost at the Rayleigh- Jeans part of the 
disc blackbody spectrum, therefore the disc surface bright- 
ness should scale as Sk oc T~ . As the outer disc temper- 
ature dominated by irradiation scales as r out oc L^/ i R~^/ 2 
(see Eq. El, having in mind R out oc a and a oc P 2 / 3 (from 
Kepler's law) as well as Lk oc SxRoutt t ne total IR lu- 
minosity of the accretion disc should be proportional to 

It is worth noting, that the above relations are simplis- 
tic, e.g. the noticeable viscous heating in the disc or con- 
tribution from the secondary star are expected to modify 
them. Therefore in order to understand more correct depen- 
dencies between the absolute NIR magnitude of LMXBs, 
their X-ray luminosity and accretion disc size (i.e. orbital 
period), we have calculated a set of LMXB models for the 
range of orbital periods of 2-750 h and X-ray luminosities 
of 1 x 10 35 -4 x 10 38 erg s" 1 . 

While varying the orbital period, we have assumed that 
the secondary star fills its Roche lobe. We assumed that 
at orbital periods less then ~6 h (M < 0.6 Mq) the sec- 
ondary star is not evolved and obey mass-radius relation 



for a main sequence star Tout et al. (19961. For the orbital 
periods larger than ~6 h (i.e. for the giant companions) 
we have fixed the companion masses at 0.6 Mq value. This 
assumption in fact does not influence our results strongly, 
and is supported by the cases when masses of the donor gi- 
ants are known, because they do not exceed this value, e.g. 
Sco X-l M 2 ~ 0.4 M ([Steeghs fc Casaresp002| >, Cyg X-2 
M 2 ~ 0.6 M ( |Orosz fc Kuulkers|1999| l~ 

After we have fixed the radius of the accretion disc (at 
the value of the tidal radius) and the shape of the disc (i.e. 
power law index n in the law H oc R n ), the main parame- 
ter, which can significantly shift the correlation up or down 
on the absolute magnitude axis, is the inclination angle of 
the binary system (this effect is illustrated in Fig. Hk. Addi- 
tional shift can be caused by different value of the fraction 
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Figure 4. Absolute NIR magnitudes of a modelled binaries ver- 
sus their orbital period and X-ray luminosity, combined into Ex 
quantity. Upper points represent absolute brightness for bina- 
ries with orbital plane inclination i = 0°, in this case Sj{ = 
(Lx/ -^Edd) 0,29 P(h) ' 92 . Lower points show absolute brightness 
for binaries with inclination i = 70°, in this case T,k = 
(Lx/iEdd) 0,28 ^(h) 93 - Dashed lines show the best linear fit to 
points. 
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Figure 5. Correlation between absolute NIR magnitude of per- 
sistent LMXBs and the combination of their X-ray luminosity 
with the orbital period, T, K = (Lx/£Edd)°- 29 -P(h) ' 92 - Crosses 
denote positions of known persistent LMXBs, dashed lines show 
best-fitting correlations calculated from our modelling, lower line 
— for inclination i = 70°, upper line - for i = 0°. Solid line shows 
the correlation Mk = 2.66 — 2.5 log 



4.2 Sample of known LMXBs 



of emission reprocessed in the accretion disc-corona system. 
We have assumed it to be 0.25 above. 

Among all models, covering the above-mentioned pa- 
rameter space we have accepted only those, which have the 
disc temperature at outer radius not lower than Th ~ 6500 K 
(see Eq. |3| . This selection was applied because we are inter- 
ested only in persistent LMXBs, while LMXBs with lower 
outer disc temperatures should be subject to disc thermal- 
viscous instability (see e.g. Lasota|2001 for review). 

The obtained absolute NIR magnitudes were fitted as 
a linear combination of logP and logLx parameters. We 
have minimized the root mean square deviations from the 
modelled Mk values. The best-fitting relation for LMXBs 
with the inclination of i — 0° is Mjf >m odel,i=o = 2.62 — 
0.731og(L x /L Ed d) - 2.29 log P(h); for the inclination of 

1 = 70°: M K ,modei,i=70 = 3.71 - 0.701og(L x /L Edd ) - 
2.32 log P(h). For both cases the rms scatter of the data 
points from the approximation is ~ 0.03. Here we adopted 
Eddington luminosity for 1.4 Mq neutron star L^dd = 

2 x 10 38 erg/s. 



Following the idea of van Paradijs & McClintock ( 1994 1 



we construct the quantity Y,k in such way that it is pro- 
portional to NIR luminosity Lk of the binary (but not 
to its absolute magnitude Mk = —2.5 log Lk + const). 
Therefore, using the average values of the obtained best- 
fitting parameters we get Ek for the near-infrared K band: 
Ek = (Lx/L E dd) ' 29 -P(h) ' 92 . For such parametrization the 
absolute brightness Mk = — 2.51ogE_R- + const. Note, that 
derived scaling is quite similar to our simplest theoretical 
estimates above. 



Having in hand physically motivated value of E/f, on 
which the LMXB luminosity should depend, we constructed 
log Ek - M K relation for the list of known LMXBs with NIR 
brightness measurements. 

We have compiled a sample of all known LMXBs with 
measured values of their brightness in K band (see Table E|. 
We have selected only binaries with neutron stars in order 
to reduce possible additional scatter due to compact object 
mass uncertainties. 

The relationship between absolute magnitude in K 
band and log Ek for known binaries from Table [4] is shown 
in Fig. [5] 

We have used linear least squares to fit this dependence 
without taking into account measurement uncertainties be- 
cause of the possible physical dispersion of the absolute NIR 
magnitude values due to the unknown inclinations of bina- 
ries. 

The best-fitting approximation of this relation for all 
considered binaries is M K = (2.78 ± 0.24) - (2.60 ± 
0.11) log Sx (1°" confidence intervals). In order to estimate 
the confidence intervals here (and below) we assumed the 
uncertainties of data points corresponding to the unity value 
of the resulted x 2 /d.o.f. Confidence intervals then were cal- 
culated as usual from interval of parameter which gives 
Ax 2 = 1. 

Note that if our parametrization of this relationship is 
correct, then the coefficient in front of logEx should be 
equal to 2.5 (because the constructed value Ek is propor- 
tional to the IR luminosity Lk of the system, whereas ab- 
solute brightness is Mk = —2.5 log Lx + const), which is 
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Table 4. Persistent LMXBs which were used to illustrate Ex — Mk correlation plotted in Fig. [5] Literature references are given where 
necessary. X-ray luminosity Lx is taken from |Revnivtsev et al.| | |2011^ who calculated it using the X-ray flux from Uhuru catalog [Forman| 
|et al.| Jl978| and the distance estimates available from the literature; orbital periods P are adopted from |Ritter fc Kolb| ( |2003[ ) ■ Values 
of Aji for 3 sources with detailed SED analysis in this paper were calculated from Ay values assuming Ax = O.llAy Values of Ak for 
remaining sources were taken from 3D Galaxy extinction map bylMarshall et al.N2006| if not indicated another reference. Magnitude of 
4U0614+091 in the K band was recalculated from J measurement of |Migliari et al.| | |2010[ | taking into account J — K colours from the 
discussed model (see the text). 



Name log L x (2-10 keV) d P m c ° TV A K 

erg s _1 kpc h Vega mag mag 



Sco X-l 


38.3 


2.8 ± 0.3 1 


18.94 


11.88 


0.1 


Cyg X-2 


38.3 


11.6 ± 0.3 2 


236.27 


12.9 


0.15 


GX 349+2 


38.2 


8.5 3 


22.5 


14. 1 4 


0.45 


GX 13+1 


37.7 


7 ± l 5 


601.7 


10. 5 6 


1.8 7 


4U1624-49 


37.5 


15.0 ± 2.9 8 


20.9 


15.9 9 


2.4 


4U1735-44 


37.7 


9.1 10 


4.65 


16.6 89 


0.16 


4U1636-53 


37.4 


5.9 


3.79 


15. 9 11 


0.27 11 


4U0614+091 


36.5 


3.2 


0.81 


16.9 


0.22 



(1) - |Bradshaw, Fomalont, fc Geldzahier]l|1999|, (2) - |Smale|||1998|, (3) - Galactic 
Center dista nce, (4) - |Wachter fc Margon ( 19 96), (5) - | Bandyopadhya y~~et al.| 
l|1999|, (6) - |Bandyopadhyay et al.| ||2002|, (7) - [ Charles fc Naylor| (|1992^ 
|Xiang, Lee, fc Nowak| (|2007| ) , (9) -|Wachter et al. " |2005| l, (10) - |Augusteijn et al.| 
l |1998| , (11) - |Russell et al.| j20li) 



compatible with the results of the fit. Fixing this coefficient 
we obtain: 



M K = (2.66 + 0.11) - 2.5 log Ex 

RMS scatter of the observed points from the fitted 
straight line for our sample of sources is ~0.3 mag. Note, 
that this value is approximately a lower limit for any Mk = 
/(Lx, P) model which does not take into account explicitly 
the binary system inclination and the time variability of its 
flux. 

As a useful application we also present here the ap- 
proximation of V — K and J — K colours, obtained with 
our model, as a function of P(h) and Lx for orbital plane 
inclination i — 0°. The best fit was calculated minimizing 
the root mean square residual on logarithmic scale: 

log(V-if+1.12) = -0.271ogLx/L E dd+0.341ogP(h)-0.75 



log(J-i i +0.28) = -0.271ogL x /L E dd + 0.331ogP(h)-1.31 

Residuals of these fits do not typically exceed 0.02-0.04 
in logarithmic scale (see Fig. Rj}. 



5 SUMMARY 

In our paper we studied relations between the X-ray lumi- 
nosity, orbital period, and absolute NIR magnitude of per- 
sistent low-mass X-ray binaries. 

We have demonstrated that LMXB spectral energy dis- 
tribution in the optical-NIR spectral range can often be ad- 
equately described by a simple model where all the surface 
elements in a binary system (comprising an accretion disc 
and a secondary star) emit black body radiation with the 
local temperature, which is itself defined by the illuminat- 
ing X-ray flux. According to this model we can construct the 



quantity Ex = (£ x /LEdd)°' 29 -P(h) ' 92 on which the infrared 
luminosity of LMXBs depends almost exclusively. Therefore 
we were finally able to define the clear relation between Ex 
and the absolute magnitude of LMXBs in the K band Mk ■ 

M K = (2.66 ± 0.11) - 2.5 log Ex 

Due to the fact that for the majority of persistent 
LMXBs we do have estimates of their X-ray luminosities, 
the essence of this relation is that it connects orbital pe- 
riod and absolute NIR magnitude of a binary system, hence 
one may consider it to be effectively the period-luminosity 
relation in the NIR. 

In some cases, however, we do see a significant emis- 
sion excess in the NIR with the respect to the prediction of 
this model suggesting the presence of an additional spectral 
component (likely optically thin synchrotron emission from 
non-thermal electrons). We therefore recommend checking 
colours in the NIR spectral bands against expected values 
using the formula we give for an adequate usage of the pre- 
sented Ex — Mk relation. Extremely red colour indexes are 
not consistent with the colour of an optically thick regions 
in the binary systems considered in this work. 

We propose that presented period-magnitude relation 
can be widely used in forthcoming surveys of the Galaxy in 
X-ray and NIR spectral domains. 
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Figure 6. Approximation of LMXB colors with models, described 
in the text. Upper plot shows V—K color, lower plot - J—K color. 
Lower panels on each plot represent residuals from the adopted 
model. 
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